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This  Interim  Technical  Progress  Report  covers  the  work  performed 
under  Contract  AF33(657)1  1265  from  January  1,  1964  to  April  1,  1964. 

It  is  published  for  technical  information  only  and  does  not  necessarily 
represent  the  recommendations,  conclusions  or  approval  of  the  Air  Force. 

This  contract  with  Republic  Aviation  Corporation  of  Farmingdale,  New 
York,  was  initiated  under  ASD  Manufacturing  Technology  Laboratory 
Project  8-1  I  I,  "Electro  Spark  Extruding.  "  It  is  administered  under  the 
direction  of  Mr.  T.  S.  Felker  of  the  Basic  Industry  Branch  MATB, 
Manufacturing  Technology  Laboratory,  Aeronautical  Systems  Division, 

W right-Patte rson  Air  Force  Base,  Ohio. 

Mr.  J.  H.  Wagner  of  the  Manufacturing  Research  Department,  Republic 
Aviation  Corporation  is  the  engineer  in  charge  of  the  project.  Mr.  Gunther 
Pfanner  is  cooperating  in  the  research. 

The  transient  wave  analysis  presented  in  this  report  was  done  by  Mr.  B.  P. 
Leftheris  of  Re-entry  Simulation  Laboratory,  Research  Division. 

The  primary  objective  of  the  Air  Force  Manufacturing  Methods  Program 
is  to  increase  producibility,  and  to  improve  the  quality  and  efficiency 
of  fabrication  of  aircraft,  missiles  and  components  thereof.  This  report 
is  being  disseminated  in  order  that  methods  and/or  equipment  developed 
may  be  used  throughout  industry,  thereby  reducing  costs  and  giving 
"MORE  AIR  FORCE  PER  DOLLAR.  " 

Your  comments  are  solicited  on  the  potential  utilisation  of  the  information 
contained  herein  as  applied  to  your  present  or  future  production  programs. 
Suggestions  concerning  additional  manufacturing  methods  development 
required  on  this  or  other  subjects  will  be  appreciated. 
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INTRODUCTION 


Electro-spark  extrusion  is  a  metalforming  concept  involving  re¬ 
covery  of  the  energy  derived  from  a  rapid  discharge  of  capacitor  stored 
electrical  power  by  reduction  of  a  metal  billet  through  an  appropriate  die. 
One  method,  the  electrohydraulic  reduction  of  the  billet,  places  the 
extrudable  material  in  an  ultra -high  pressure  fluid  environment  with  the 
rapid  energy  lischarge  occurring  across  a  suitable  electrode  gap.  This 
short  time  duration  gap  discharge  in  a  properly  designed  high  strength 
container  causes  a  shock-pressure  wave  to  strike  the  billet  face  and, 
upon  reflection,  to  initiate  partible  acceleration  in  the  billet  material 
and  to  release  to  the  billet  face  some  of  the  intrinsic  kinetic  and  pressure 
energy  available  as  a  consequence  of  known  wave  reflection  parameters. 
The  sum  of  these  phenomena  when  added  to  the  potential  energy  present 
in  the  high  fluid  pressure  environment  have  been  shown  by  mathematical 
analysis  to  be  ample  under  carefully  controlled  conditions  to  cause  a 
metal  billet  to  yield  and  to  extrude  a  given  length  for  the  time  duration 
of  one  electrical  pulse.  The  problem  of  maintaining  extrusion  then 
properly  becomes  one  of  providing  sufficient  discrete  pulse  discharges 
with  the  correct  time  interval  to  allow  a  constant  velocity  extrusion. 

In  addition,  the  capacitor  discharge  equipment  can  be  alternately 
employed  to  produce  impulsive  electromagnetic  forces  upon  the  billet  in 
the  hydrostatic  container.  In  this  method,  the  rapid  discharge  of  the 
sfc~cd  electrical  energy  produced  high  magnetic  fields  about  an  appro¬ 
priately  constructed  conductive  coil.  By  induction,  to  a  conductive 
pistoi. -disc  interposed  between  the  coil  and  the  billet,  lepulsive  magnetic 
forces  rapidly  accelerate  the  disc  and  thereby  drive  a  pressure  wa  : 
against  the  billet. 

As  an  alternative,  where  the  billet  material  is  sufficiently  conductive 
the  electromagnetic  repulsion  may  be  coupled  to  the  billet  face  directly  in 
order  to  induce  the  acceleration  forces  required  to  extrude  the  billet. 

A  third  method  of  capacitor  stored  energy  utilization  is  to  recover 
the  energy  created  by  the  electrohydraulic  discharge  and  to  transfer  it 
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by  mechanical  means  to  the  billet.  Displacement  of  the  billet  is  accom¬ 
plished  mechanically  by  the  use  of  a  piston-ram  on  the  billet  material  in 
a  suitable  container.  The  pressure  pulse  or  electromagnetic  energy 
created  by  a  series  of  discharges  impinges  on  the  face  of  the  piston  causing 
it  to  move  forward  to  extrude  the  billet. 

The  advantages  offered  by  the  foregoing  approach  to  electrohydraulic 
extrusion  are:  (1)  elimination  of  billet  container  wall  friction,  (2)  re¬ 
duction  of  billet  die  surface  friction,  (3)  creation  of  pressures  consider¬ 
ably  beyond  those  obtained  hydrostatically,  (4)  reduced  container  size  and 
cost. 


To  determine  the  production  potential  for  extruding  a  steel  alloy  by 
the  use  of  capacitor  stored  energy,  the  Aeronautical  Systems  Division 
of  Wright  Patterson  Air  Force  Base  has  awarded  Contract  No.  AF33(656)- 
11265  to  Republic  Aviation  Corporation.  Extrusion  will  be  attempted  in 
the  course  of  the  program  by  conducting  experiments  with  capacitor  dis¬ 
charge  equipment  to  develop  the  most  suitable  techniques  for  metal  reduc¬ 
tion  by  several  electrohydraulic  methods.  The  program  consists  of  two 
phases  as  follows: 

Phase  I  -  Design  of  the  Extrusion  Equipment 
Phase  II  -  Development  of  the  Extrusion  Process 

The  approaches  mentioned  in  the  foregoing  are  to  be  Thoroughly 
investigated  for  suitability  fward  attaining  the  objectives  of  the  program; 
i.  e. ,  extrusion  of  a  2  inch  diameter  to  a  1/4  inch  round  . 

Mathematical  treatment  of  the  electrohydraulic  pulse  extrusion  case 
(Reference  1)  produced  a  descriptive  equation  containing  a  combination  of 
dynamic  and  static  variables  generally  similar  to  the  well  known  equation 
used  for  static  extrusion  (  P  =  K  In  -j- ).  In  the  static  extrusion  case, 
the  pressure  is  raised  until  extrusion  starts;  while  in  the  pulse  extrusion 
case,  the  pulse  propagates  through  the  billet  imparting  a  motion  to  the 
billet.  The  extrusion  produced  for  one  pulse  depends  upon  the  pressure 
amplitude,  while  the  final  extruded  length  depends  upon  the  rate  the 
pulses  are  generated. 
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From  the  foregoing,  it  is  apparent  that  achievement  of  the  total 
energy  necessary  to  extrude  a  target  section  requires  multiple  electrical 
discharges.  The  time  between  discharges  should  be  minimum  to  reduce 
the  overall  extrusion  time  and  avoid  stop-start  extrusion.  In  order  to 
optimize  these  requirements,  it  is  necessary  to  modify  the  156, 000  joule 
Republic  Aviation  capacitor  discharge  facility  to  include  the  necessary 
equipment  to  allow  faster  charging  rates  and  to  design  a  switching 
arrangement  and  duty  cycle  timer  to  favor  the  rapid  cyclical  release  of 
stored  energy  at  precisely  timed  intervals. 

The  theory  of  unsteady  waves,  (Reference  1)  used  to  develop  working 
equations  for  computing  the  radial  and  hoop  strains  in  hollow  cylinders 
under  transient  internal  loads  indicated  that  pressure  vessel  design 
requirements  to  contain  dynamic  forces  superposed  over  high  hydrostatic 
pressures  was  possible  and  practicable.  This  analysis  pointed  up  desirable 
container  dimensional  geometry  and  indicated  an  order  of  magnitude  of 
mechanical  properties  for  the  materials  in  the  critical  areas  of  the  chamber. 
These  criteria  have  been  formulated  into  general  design  requirements  and 
are  presented  in  this  report. 
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DESIGN  OF  A  PRESSURE  VESSEL 
FOR  ELECTROHYDRAULIC  EXTRUSION 


A.  Basic  Equipment  Design 

1.  The  Spark  Discharge  Extrusion  Concept 

To  accomplish  extrusion,  power  is  initially  supplied  to  a  capa¬ 
citor  bank  and  the  stored  charge  held  off  by  a  suitable  switch  until  re¬ 
leased  to  the  high  pressure  device  containing  an  open  gap  electrode,  an 
extrusion  die  and  a  metal  billet.  The  billet  is  maintained  under  high 
hydrostatic  pressure  defined  by  the  mechanical  properties  of  the  part¬ 
icular  extrudable  material  used  and  the  degree  of  metal  reduction  required. 
Once  the  storage  of  energy  is  completed,  the  discharge  event  begins  when 
a  trigger  switch  {vacuum  gap)  is  ionized  and  current  flows  into  the  circuit, 
thereby  ionizing  the  high  fluid  pressure  gap  between  the  electrodes  (refer 
to  Figure  1).  A  spark  channel  is  abruptly  created  and  the  vapor  products 
in  the  channel  expand  into  a  spherical  "gas  bubble. "  The  inertia  of  the 
circumjacent  fluid  coupled  with  the  extremely  high  velocity  of  expansion 
of  the  "gas  bubble"  produces  a  compression  shock  wave  in  the  surrounding 
liquid.  This  compressed  water  layer  then  travels  through  the  liquid  as 
a  shock  pressure  front  at  approximately  acoustic  velocity  with  a  time 
duration  of  approximately  20  microseconds. 

The  energy  relea  :d  by  the  gap  discharge  consists  of  kinetic  energy 
plus  the  energy  due  to  the  increased  pressure  in  the  traveling  wave.  At 
the  billet  face  the  pressure  approximately  doubles  because  of  the  rein¬ 
forcement  of  the  incident  and  reflected  waves.  The  shock  pressure  wave 
splits  into  a  reflected  pulse  back  into  the  fluid  and  a  transmitted  pulse 
forward  into  the  billet  material.  The  billet  attains,  by  the  same  phenomena 
as  the  fluid,  a  particle  velocity  sufficient  to  cause  longitudinal  movement 
for  a  length  proportional  to  the  energy  released. 

The  very  benefits  derived  from  the  shock  pressure  pulse  in  extrusion 
(i.  e. ,  particle  velocity  sufficient  to  cause  plastic  strain  and  movement  of 
the  billet)become  a  limiting  factor  when  applied  directly  to  the  structure 
in  the  pressure  vessel  walls.  As  the  cylinder  is  initially  under  high 
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internal  static  pressures  the  effects  of  the  superposed  dynamic  loading 
become  critical.  Results  of  the  analysis  completed  in  subject  Reference  1, 
indicated  a  practical  container  geometry  is  attainable  to  withstand  the 
explosive  forces  from  an  electrical  gap  discharge  event. 

The  illustration  (Figure  1)  shows  the  spark  discharge  extrusion 
process  and  equipment  in  a  conceptual  manner.  A  capacitor  bank  to 
supply  the  energy,  a  high  pressure  chamber,  pumping  system  and  asso¬ 
ciated  controls  are  shown. 

2.  Effect  of  Shock  Loading  in  Container  Cavity  Dimension 

Equations  have  been  previously  derived  and  calculations  prepared 
in  Reference  1  for  the  general  problem  of  a  thick-walled  cylinder  under  the 
influence  of  a  transient  internal  load  where  the  unbalanced  forces  from  a 
detonative  capacitor  discharge  set  up  velocity  gradients  in  part  of  the 
cylinder  while  the  remainder  is  completely  at  rest.  The  problem  has 
been  formulated  and  explicit  solutions  given  for  elastic  and  elasto-plastic 
cases.  These  specialized  interpretations  are  presented  in  Appendix  I 
for  reference. 

With  the  aforementioned  analysis  in  mind,  a  tentative  pressure 
vessel  construction  (Figure  2)  using  a  single  shrink  shroud  of  4340  steel 
(24  inch  O.  D. )  over  an  18%  Ni,  maraging  steel  liner  (2  inch  I.  D. )  with  a  static 
internal  pressure  of  200,000  psi  was  derived  using  the  customary  Lame 
equations.  A  desirable  residual  compressive  hoop  stress  resulting  from 
the  combination  of  shrink  .it  assembly  and  cold  work  fabrication  proced¬ 
ures  in  the  most  highly  stressed  portion  of  the  apparatus,  (tV>  inside 
diameter;  rQ  =  1  inch)  was  required  to  withstand  safely  hydrostatic  internal 
pressures  io  the  maximum  amount  of  200,000  psi. 

The  question  of  a  serviceable  design  to  contain  the  forces  re¬ 
leased  by  the  discharge  of  electrical  capacitor  stored  energy  in  addition 
to  the  loading  imposed  on  the  equipment  by  the  high  order  of  hydrostatic 
oressure  presented  a  problem  to  which  no  exact  solution  could  be  found. 

For  this  reason,  the  mathematical  analysis,  Reference  (l),was  prepared 
to  predict  the  effects  of  detonative  shock  loads  on  a  container  using  normal 
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construction  materials  of  significant  mechanical  strength  values. 

This  analysis  indicated  that  residual  hoop  stresses  result  in 
the  container  walls  from  reaction  of  the  metal  to  the  passing  of  a  maximum 
detonative  shock  pressure  pulse.  The  maximum  response  occurs  at  the 
fluid  metal  interface  at  the  internal  radius  (r0  =  1  inch)  where  the  pressure 
pulse  is  greatest.  The  material  has  been  shown  by  the  analysis  to  exper¬ 
ience  first  plastic  flow  in  the  outward  direction  and  then  inward,  thus 
giving  rise  to  the  possibility  of  secondary  yielding. 

Reference  to  Case  1  in  Aopendix  I,  the  hoop  strain  for  zero 
stress  (permanent  set)  is  0.  0036,  The  residual  strain  after  unloading 
is  0.  00463  (compression).  The  diagram  for  this  set  of  conditions  indicates 
that  a  permanent  set  to  a  residual  compressive  stress  level  at  the  bore 


ol'  453,  000  psi.  This  is  a  value  greater  than  the  ultimate  strength  of  the 
liner  material  (250,000  psi)  being  considered  without  the  addition  of 
hydrostatic  pressure.  Therefore,  the  assembly  shown  in  Figure  2  is 
proven  unsatisfactory,  because  of  the  high  degree  of  deformation  present 
et  the  critical  location  vQ  =  1  inch  from  the  action  of  the  passing  of  the 
shock  wave. 

As  the  reverse  yielding  at  rQ  =  1  inch  has  been  shown  to  be  beyond 
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the  mechanical  property  values  of  18%  Ni.  maraging  steels,  it  is  of  interest 
to  find  a  radius  where  the  severity  of  the  shock  wave  has  been  sufficiently 
diminished  so  as  to  present  a  less  drastic  effect  on  the  walls  of  the  pressure 
container.  A  plot,  Figure  3  of  the  particle  velocity  against  the  radius  in 
the  container  wall  indicates  that  the  particle  velocity  attenuates  rapidly  from 
its  initial  critical  value  of  u  =  2864  inches  per  second.  It  is  reasonable  to 
assume  that  a  wall  radius  can  be  found  where  the  effect  of  the  particle 
velocity  on  the  container  walls  is  reduced  to  the  degree  that  severe  deform¬ 
ations  do  not  occur  from  the  loading-unloading  cycle  of  the  transient  wave. 
The  curve  changes  slope  between  the  values  r  =  2  and  r  =  3  inches.  At 
r  =  3  inches,  a  40%  attenuation  has  diminished  the  particle  velocity  to  a 
value  of  1654  inches  per  second.  This  value  of  rQ  =  3  inches  was  selected 
for  analysis. 

The  calculations  for  an  rD  =  3  inches  are  outlined  in  Appendix 
I,  Case  (2). 

The  hoop  strain  imposed  at  r0  =  3  inches  from  the  action  of  the 
loading  wave  where  the  pressure  has  attenuated  to  248,280  psi  has  been 
shown  using  the  concept  of  elastic  recovery  to  be  equal  to  zero,  thus 
indicating  no  permanent  set  in  the  hoop  direction  from  the  loading  wave 
action.  Due  to  the  unloading  there  is  elastic  recovery  strain  in  the  amount 
of  0.  00368  (compression).  The  equivalent  stress  for  this  strain  value  is 
110,000  psi  (compression).  The  significance  of  this  is  that  a  permanent 
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residual  hoop  stress  (compression)  due  to  the  action  of  the  wave  remains 
at  the  inside  diameter  of  the  chamber  (rQ  =  3  inches).  If  there  were  no 
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considerations,  other  than  hoop  stress,  the  pressure  vessel  at  rQ  =  3  inches 
could  elastically  withstand  dynamic  stresses  from  -110,000  psi  to  the  yield 
point  of  the  material,  250,000  psi  or  a  total  of  360,000  psi  in  the  hoop 
direction. 


It  must  be  remembered  that  the  analysis  outlined  in  the  Appendix 
I  indicates  the  condition  in  the  pressure  vessel  walls  at  rQ  =  3  inches.  This 
case  is  for  the  same  wave  (at  a  later  period  in  time)  as  was  analyzed  for 
rQ  =  1  inch.  The  particle  velocity  and  the  pressure,  of  course,  are  atten¬ 
uated  from'passage  through  the  material  of  the  pressure  vessel  walls. 
Although  the  calculations  show  a  small  residual  compressive  stress  of 
30,000  psi,  this  value  can  be  disregarded  if  we  consider  rn  =  3  inches  to 
be  the  container  cavity  dimension  at  the  fluid-metal  interface.  (This  will 
be  made  further  evident  in  the  following  section.! 

3.  Combined  Effect  of  Internal  Detonative  and  Hydrostatic  Loading 

Because  the  shock  wave  loading  by  itself  has  been  previously 
shown  to  behave  elastically  between  1 10,  000  and  the  yield  point  of  the 
material  in  the  wall  at  rQ  =  3  inches,  it  becomes  of  interest  to  subject 
the  stresses  produced  from  a  combination  of  the  dynamic  and  static 
loading  to  analysis  to  see  if  the  coupling  can  be  accomplished  without 
failure  in  the  wall.  The  calculations  for  this  are  outlined  in  Appendix  II. 

The  work  indicates  that  a  compressive  pre-stress  in  the  hoop 
direction  of  approximately  133,000  psi  must  be  incorporated  in  the  pressure 
vessel  construction  by  some  method  to  insure  that  failure  will  not  occur 
at  rQ  =  3  inches  as  a  result  of  the  combining  of  the  hydrostatic  and  dynamic 
loading.  This  order  of  compression  stress  is  to  be  obtained  by  the  use  of 
careful  shrink  fit  procedures  in  addition  to  appropriate  autofrettage  of  the 
inner  liner  prior  to  assembly.  As  an  interference  shrink  fit  construction 
will  be  used.  The  interference  must  be  calculated  commensurate  with 
ordinary  assembly  procedures.  In  addition,  the  residual  stress  level 
after  shrink  and  at  pressure  must  be  maintained  so  that  this  area  jvill  not 
be  subject  to  failure  during  usage. 
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4.  Replaceable  Liner  Construction 


The  analyses  and  calculations  presented  for  r0  =  1  inch  and 
rQ  =  3  inches  container  dimensions  have  considered  the  decay  of  one 
transient  detonative  shock  wave  in  a  cylinder  from  a  point  at  the  pressur¬ 
ized  fluid-metal  interface  outward  for  a  distance  of  three  inches.  Figure  4 
shows  the  respective  condition  of  one  wave  at  rQ  =  1  inch  and  rQ  =  3  inches 
at  two  different  instants  of  time.  If  we  should  consider  a  chamber  only  of 
cavity  diameter  of  6  inches  (r0  =  3  inches)  the  pressure  upon  reflection 
would  again  be  430,  000  psi  and  the  effect  on  the  container  wall  would  be 
too  severe  to  contemplate  for  tho  proposed  experiments.  If  wo,  however, 
take  the  maximum  pressure  {430,000  psi)  at  the  2  inch  diameter  and  allow 
the  effects  (particle  velocity  and  pressure)  to  attenuate  by  some  means, 
the  consequences  of  the  event  at  the  6  inch  diameter  interface  can  be 
tolerated. 


The  means  proposed  to  mitigate  the  severity  of  the  shock  pulse 
on  the  wall  material  is  the  insertion  of  a  removable  liner  2  inch  I.  D.  x 
6  inch  O.  D.  in  the  pressurized  container.  This  liner  will  be  completely 
restrained  by  the  hydrostatically  pressurized  fluid  in  the  radial,  tangential 
and  axial  direction.  Such  a  construction  is  shown  in  Figure  5. 

It  is  anticipated  that  over  a  period  of  time  after  a  specific 
number  of  maximum  duty  cycles,  the  replaceable  liner  should  fall.  When 
this  occurs,  the  liner  c-  •>  be  removed  and  a  replacement  provided. 

The  initial  static  and  dynamic  loading  cycle  is  depicted  graphically 
in  Figure  6.  ^  ’s  the  initial  prestress  at  F  =  0.  Upon  the  insertion  of 
pressure  into  the  cavity,  the  stress  is  raised  to  Oi  +  (%.  The  first  maximum 
discharge  resets  this  stress  pattern  to  the  value  at  (%and  upon  release  of 
pressure,  the  residual  hoop  stress  returns  to  the  point  of  (%  This  is  then 
the  permanent  working  stress  in  the  cavity  of  the  pressure  chamber. 
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STATUS  OF  THE  WORK 


The  analyses  derived  for  the  effect  of  detonative  shock  loading  in  a 
thick  walled  cylinder  hydrostatically  loaded  have  been  completed  and  the 
conclusions  derived  therefrom  have  been  applied  to  the  design  of  a  suitable 
piece  of  equipment  to  be  used  in  the  extrusion  experiments. 

All  vendor  proposals  for  equipment  have  been  reviewed  for  conform¬ 
ance  to  specification  and  suitability  of  costs,  and  based  upon  this  review  , 
an  agreement  has  been  entered  into  with  Harwood  Engineering  Company, 
Walpole,  Massachusetts  to  develop  with  Republic  Aviation  Corporation, 
and  to  manufacture  the  equipment  required  for  the  proposed  experiments. 
”arnouu  has  started  fabrication  of  a  200,000  psi  high  pressure  pumping 
system  to  displace  20  in  3/min.  Engineering  has  been  completed  in  the 
design  of  the  pressure  vessel  and  work  has  commenced  with  the  assembly 
drawings.  These  should  be  ready  for  detailing  early  in  the  next  reporting 
period. 

Work  has  proceeded  with  the  electrical  considerations  necessary  to 
alter  the  Republic  Aviation  Corporation  156,000  joule  capacitor  bank  to 
attain  suitable  charging  rates  and  to  release  energy  to  the  extrusion 
chamber  at  the  repetition  rates  that  are  anticipated  to  be  necessary  to 
accomplish  extrusion.  This  work  will  continue  for  the  next  reporting 
period. 
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APPENDIX  I 


CALCULATION  FOR  STRAIN  DEFORMATION 
DUE  TO  TRANSIENT  WAVE 


I  Statement  of  the  Problem 

Consider  a  thick-walled  cylinder  with  =  250,000  psi  and  the  following 
dimensions  ( the  cylinder  is  filled  with  water): 


INCHES 

FIGURE  1 

i  ;riangular  wave  of  the  following  form  is  transmitted  along  the  cylinder. 


This  wave  reflects  upon  arrival  on  a  closed  end.  Find  the  resulting  strains 

r  =  1  inch  and  r  =  3  inche. . 
o  o 

The  pressure  of  the  reflected  wave  front  is  approximately  twice  the  pressure 
of  the  incident  wave  front.  This  is  considered  the  highest  possible  internal  p>  .in¬ 
sure;  thus,  the  calculations  will  be  performed  for  the  reflected  region  only. 

To  illustrate  the  wave  decay,  two  cases  are  analyzed;  Case  1,  where 
rQ=  1  inch,  and  Case  2,  where  rQ  =  3  inches.  Hence  P=- 2 (215, 000)  =-430,000 
ost  at  rQ=  1  inch.  Using  equation  (23),  we  have 


i 
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where 


u ,  =  -  S-  =  -  2864.6  in./sec 

M>  0 

p0  =  0.29lb/in.3 
CQ  =  200,000  in./sec 
g  =  32.2ft/8ec2 

From  equation  (34),  the  velocity  at  rQ  =  3  inches  (chosen  arbitrarily)  is 
u„  =  u,  ,/~  =  -  2,864.6  JJ~=  -  1654.0  in./sec 

o  l  V  Tg  Vo 

Thus  the  pressure  at  rQ  -  3  inches  is  given  by 
P3=  JY  *  I3!  =  -  248,280  psi 


A.  CASE  1 


1. 


Loading  Wave 

From  equation  (28)  the  density  ratio  at  t  =  0  is  given  by 


2.00x  105 
2.00X105  -  2864.6 


=  .98568 


The  maximum  cTj  at  the  start  of  plastic  deformation  is  given  by 

ffy  250,000  „„„  , 

°1=  TTH  =  ‘  -o^7~=-  357,143  psi 


where 


av  =  -  250,000  psi 


l-H  =  0.7 


Thus  C=l“  (i-2x0.32)=  1-.0119X. 82  =  .99042 

30x10° 


at  ay-~  357* 143  P£i 

As  soon  as  the  material  begins  stretching  out,  the  loading  becomes  plastic. 
Considering  1  =  0  from  equation  (38),  we  have 
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«  =  -/5y.C  =  /  00833  x.  98568  =  7. 0082106  =  .09061 

where  5  =  =  .00833 

y  30 x 10® 


Also, 


2864.6 

2xl05 


-.01432 


Substituting  in  equation  (39),  we  have 

/  =  cosh  .  09061 1  -  fljUff  sinh  •0f,08lJ  ‘  1 

Since  u  =  ,  we  have 

u  =  .09061  (sinh  .  09061 t-  .15800  cosh  .  09061 1} 


First,  the  time  tQ,  where  0  =  0,  is  found  by  plotting  u  versus  t. 

(''ie  Figure  2).  The  £  versus  t  plot  is  shown  in  Figure  3:  it  is  shown  that  as 

u  approaches  zero,  £  approaches  a  constant  value  of  0.0125.  The  radial  strain 
r 

is  given  by  <j  =  C  —  "  1:  The  ci  versus  t  plot  is  Bhown  in  Figure  4:  e1  ap¬ 
proaches  the  value  of  .  0264. 


2.  Unloading  Wave 

The  duration  of  the  triangular  wave  is  20  fi  sec:  its  equivalent 
rectangular  profile  will  have  a  duration  of  10  Jjsec.  Hence, 

r  C  1  2.0xl05xl0x  10"6  „  „ 

1  ~  *r' '  "  I 

o 

It  is  shown  in  Figure  2,  however,  that  the  loading  process  terminates  at  t  =  1.75, 
Thus,  when  the  unloading  wave  occurs,  the  material  is  in  equilibrium. 

For  the  unloading  wave,  we  have  the  following  constants  (equation  (53)): 
I  =  1 

C  =  1:  The  material  returns  to  its  initial  volume. 
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Hence  a  =  \  +  M  =  °*  5  +  0.303  =  .803 

where  (i=  0.303.  N.B.  The  value  of  after  unloading  1b  ze.o 

0  =  l-n  =  0.697 

y  =  0.303 

u  =.  01432 
o 

£  =  -  1.0000  -  =  -  1.0000  -  .01432  =  -  1.01432 

o 


ThUS  S3'43472  mA  2  =°-4015 

Hence  4  +  f  (jjf)  =  “  1.01432  +  0.4015  ( .43472) 

=  -  1.01432  +.17454  =  -  .83978 

^  Y=  --S'-- W8195 

2 

Also,  x  |  =  .8585  x  .43472  =  .37321 

. 43472  x. 5  =.21736 

®  Y  =  -  .4015  x  .978195  =  -  .392745 

|  X  2  =  ,43472  x  .4015  =  .17454 

£  =  -  .978195  =  -  .489097 

i  5 

and  |  jj  -  1  =  -  1  +  *21736  =  “  .78264 
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Hence  u  =  .  2',  686  sinh  .  8585 1  -  .  79696  cosh  .  8585 1  +  .  78264  e  * ' 4015 1 

and  /  =  e",4015t  {.43472  cosh  .8585t  -  1.01432  sinh  .8565  tj-.  43472 

In  the  plot  of  u  versus  t  ,  Figure  5.  the  time  when  the  velocity  becomes  zero 
is  shown  as  tQ  =  0.026  (t=  ,26x  10~6  sec).  Substituting  in  the  equation  of  g  we 
find  the  residual  hoop  strain,  £  =  .01713,  Equations  (51)  and  (52)  were  derived, 
however,  with  g  -  0  at  t  =  0  (with  reference  to  the  equilibrium  conditions  after 
the  loading). 


To  find  the  residual  hoop  strain  after  unloading,  we  have: 

-  ^unloading  "  ^"loading 

=  .01713  -  .0125  =  jyMgg  (compression) 


In  order  to  find  the  residual  hoop  stress  we  use  the  principle  of 
elastic  recovery.  Thus,  the  volumetric  strain  ey  is  given  by: 

ev  =  *1  +  ^unloading 


where  Cj  is  the  radial  linear  strain  change  upon  relief  and  is  the  corresponding 

hoop  strain.  We  know,  however,  the  radial  unloading  stress:  i.e.,  Ac.  =  430,000 

psi.  Hence,  e.  =  430>  =  .0143. 

1  30x10° 


The  value  of  ev  is  also  known  from  the  loading  conditions: 


Thus,  the  hoop  strain  for  zero  stress  (permanent  set)  is  given  by 
fg 0  =  ev~  (l  =  • 01813  -  • 0143  =  ’0Q38‘ 

The  corresponding  permanent  set  in  the  radial  direction  is 
c  = .  026  -  .  0143  =  .012 .  In  the  calculations  for  the  unloading  wave,  it  was 
assumed  that  upon  relief  the  recovery  was  instantaneous:  thus  —  =  1  at  t=  0. 
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FIGURE  5 


Then,  for  t  >  0,  -  decreases,  changing  both  hoop  and  radial  strains  according  to 
to  the  following  equations: 


Without  pursuing  this  case  any  further,  it  is  easily  concluded  that  secondary 
yielding  will  take  place.  This  is  undesirable  for  most  designs. 

B.  CASE  2 

1  •  .J  oading  Wavp 

The  velocity  at  rQ  =  3  inches  was  given  as: 


.  %  _  .00827  =  -  .0909885 

(i)  .0908906 

Hence  /  =  cosh  .  0908906 1  -  .0909885  sinh  . 0908906 t  -  1 
and  u  -■  .0908906  {sinh  .  0908906 1  -  .0909885  cosh  .0908906  } 

The  duration  of  the  wave*  on  the  other  hand,  is  given  by 

r  10  x  10~6  x  2  x  105  2 

3  3 

where 

t  =  10  x  10"6 

CQ  =  2  x  106 

r  =3  Inches 
o 

The  u  versus  t  and  £  versus  t  plots  are  shown  in  Figures  6  and  7, 
respectively.  Another  case  of  cry  =  150, 000  psi  is  also  shown  for  comparison. 

With  the  value  of  t  at  rQ  =  3  inches  being  smaller  than  at  rQ  =  1  inch, 
it  is  shown  that  the  unloading  wave  occurs  before  the  material  achieves  equilibrium. 

2.  Unloading  Wave 

Following  the  same  procedure  as  for  rQ  =  3  inches,  with 

u„  =  0.00827  -  .00275*  =  .00552 
o 

$  =  e",4016t  {.43472  cosh  .8585t  -  .96794  sinh  ,3585t  }  -  .43472 

and  fl  -  ~  .  27793  sinh  .  8585 1  +  .  78816  cosh  .  8585 1  -  .  78264  e* 401 5  * 

The  plot  of  Q  versus  t  is  shown  in  Figure  8.  Comparing  this  un¬ 
loading  with  that  of  rQ  =  1  inch,  we  see  that,  at  rQ  =  3  inches,  equilibrium  is 
achieved  in  about  half  the  time. 


*  velocity  from  the  loading  process 
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The  corresponding  value  of  £  at  t  =  .01  is  £  =  .  00727.  The  <f  value 
at  t  =  0.667  for  loading  is  £  =  .00365.  Hence,  the  resulting  strain  is  (.00733 
-  .00365)  or  ^  =  ,00368  (compression). 

To  find  the  resulting  stresses,  the  permanent  set  in  both  the  hoop 
and  radial  directions  must  be  found.  Employing  the  concept  of  elastic  recovery 
again,  we  have 

.  ev=  ^>  +  <l 

f  .  248.2M  .  .008276 
30x10° 

Furthermore,  -e=  +  -^=-fi_+i  =  -.  99173  +  1  =.  00827. 

O  p0 

Thus,  Q  =  0:  This  result  indicates  that  there  is  no  permanent  set  in  the  hoop 

direction.  In  the  radial  direction,  ~  =  l  +  f  =  1+  .00365  =  1.00365  from 

o 

loading.  Hence  ~  =.99636.  Also,  ^  =  .99173. 

*1  =  ‘  99036  x  .  99173  -  1  =  .  98812  -  1  =  -  .  01188  from  loading. 

During  the  unloading,  there  is  the  elastic  recovery  where  e^  =  .  00827: 
then,  there  is  the  further  change  of  Cj  due  to  the  —  change. 

1  ro 

r 

Hr;uce  =~  ■  It-  .00827 

§"  =1-^-1-  .00727  =  .99373 
"o 

r 

~  =  1.0063 

and  (x  =  .0063  +  .00827  =  .01457 
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Thus,  change  in  =  .01457  -  .0119  -  .00265  (tension).  The  value 
of  when  the  radial  stress  is  zero  is  -.01190  +  .00827  -  -.0036. 

Residual  radial  stress  =  -  .0036  +  .0026  =  -  .001  or  -  .001  x  30  x  106 
=  -30, 000  psi  Compression). 

The  residual  hoop  stress  =  .00368  x  30  x  106  =  110,400  psi 
(compression).  From  the  secondary  yielding  criterion,  we  have 

cr2  -  Oi  <  ay 

or  -  30, 000  +  no,  000  =  80, 000  psi.  Since  ay  was  assumed  at  250, 000  psi,  the 
material  will  not  go  through  secondary  yielding. 


APPENDIX  II 


CALCULATIONS  SHOWING  COMBINED  EFFECT 
OF  IMPULSE  AND  STATIC  LOADING 


I.  Introduction 

The  combined  effect  of  the  action  of  a  transient  shock  wave  over  an 
internal  hydrostatic  pressure  in  thick  walled  cylinder  is  analysed  in  the 
following  to  verify  the  structural  integrity  of  a  container.  The  method  of 
superposition  of  stress  is  used. 

If  we  are  given  a  cylinder  with  the  following  dimensions: 


a  =  3" 
b  =  15" 
w  =  b/a  =  5" 

P  =  internal  pressure 
(200,  000  psi 


For  the  purpose  of  discussion,  we  define  the  following  variables  at 
an  internal  radius  of  r0  =  a  =  3  inches: 

<Ti  =  hoop  stress  under  internal  hydrostatic 
pressure  alone 

c%  =  radial  stress  under  internal  hydrostatic 
pressure  alone 

or3  =  residual  hoop  stress  in  bore  at  zero  pressure 
0,  =  residual  hoop  stress  resulting  from  an  impulse 

C%  =  03  +  04  =  resultant  hoop  stress  after  release  of 
internal  pressure  due  to  passage  of 
first  maximum  shock  impulse 
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c%  =  ff\  +  +  <r4  =  combined  hoop  stresses  under 

static  pressure  after  passage  of 
impulses 

II  Calculation  for  Stress  from  Internal  Pressure  Alone 

The  equations  of  Lame  may  be  used  to  calculate  the  elastic  stresses 
in  a  thick  walled  cylinder,  internally  loaded.  The  relationship  for  the 
hoop  stress  at  the  bore  states  that. 


Then,  for  a  wall  thickne s  w  =  5  and  an  internal  pressure  of 
P  =  200,  000  psi 

Oi  =  200,000  =  217,000  psi  (tension) 

By  definition, 

o^=  -P  =  -  200,000  psi  (compression) 

From  the  secondary  yielding  criterion,  we  have 

oy  >  os  - 

or  -200,000  -  217,000  =  -417,000  psi.  Since  oy  was  assumed  at 
250,000  psi,  the  material  will  go  through  secondary  yielding. 

However,  to  reduce  the  stress  to  a  point  under  the  250,000  psi 
yield  stress  of  the  material,  we  must  impose  a  negative  hoop  stress  in 
the  bore  prior  to  loading.  To  bring  the  bore  to  a  zero  hoop  stress  at  an 
ji/terr  pressure  of  200,000  psi  let, 

<r3  =  -2 17,  000  psi 

Then  using  the  Tresca  Criterion,  we  have 
oy  >  Os  -Oi 

or  -200j  000  +  0  =  -200, 000  psi.  This  value  is  suitably  below  the 
required  250,000  yield  stress.  Therefore,  the  material  will  not  fail  and 
we  can  superpose  the  dynamic  loading  to  predict  the  final  state. 

Ill  Continued  Dynamic  and  Static  Load 

Previous  analysis  (Appendix  I)  has  shown  that  the  residual  hoop 
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stress  a,  after  discharge  (peak  pressure,  P  =  248,280  psi  and  time 
duration  t  =  20  /i sec  is, 

c r4  =  -110,000  psi 

The  combined  hoop  stresses  (  ct6)  under  static  pressure  after  the 
initial  capacitor  discharge  will  then  be 

(%=  hi  +<r3)  + 

•  =  217,000  -217,000  +  (-110,000) 

=  -110,000  psi  (compression) 
the  radial  stress, 

=  -P  =  200,000  psi  (compression) 

From  Tresca's  Criterion,  we  have, 

Oy  >  oa  -  Ol 

or  -200,000  -  (-110,000)  =  -90,000  psi.  Since  ay  is  again  250,000  psi 
yielding  will  not  occur  in  the  bore  from  the  combined  stresses  due  to 
dynamic  and  static  forces. 

IV  Release  of  Pressure 

Upon  the  release  of  the  static  pressure,  the  permanent  hoop  stress 
imposed  by  the  first  discharge  event  will  reverse  by  the  amount  of  Oi. 

Then  the  ~esidual  stress  with  internal  pressure  P  =  0  will  be, 

Ob  =  Oa  +  O4 

=  -217,00'J  -110,000 
=  -317, 000  psi  (compression) 

Since  cg  exceeds  the  250,000  psi  yield  strength  of  the  material 
under  the  condition  of  no  internal  pressure,  the  bore  will  yield  radially 
inward.  This  is  an  unsatisfactory  stress  condition  for  cylinder  design. 

V  Solution  of  the  Problem 

Because  yielding  has  occurred  at  the  bore  with  an  initial  hoop  stress, 
cr3  =  217,000  psi  (compression)  a  03  must  be  found  that  will  allow  elastic 
behavior  of  the  material  within  the  250,000  psi  yield  strength  wall  material 
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being  considered* 

Let  us  calculate  a  new  residual  hoop  stress  in  the  bore  at  no 
pressure  that  will  allow  the  desired  elastic  action  within  the  limits  of 
the  250,000  psi  yield  strength  material.  Then 

ofe  =  at  +  Ob  + 

240,000  =  217,000  +  <%  +  (-1 10,  000 ) 

Ob  =  -133,  000  psi. 

Therefore,  the  desired  value  of  pressures  at  the  bore  =  133,  000  psi 
(compression) 

The  combined  stress  (r%)  at  a  pressure  P  =  200,000  after  a  maximum 
energy  impulse  is  then 

=  Oi  +  o  3  +  ^4 

=  217,000  +  (-133,000)  +  (-110,000) 

=  -26,  000  psi  (compression) 
oj,  =  -P  =  -200,  OOOpsi  (compression) 

Then  using  the  Tresca  Criterion, 
cry  >  <ra  -  ffj 

or  -200,000  -(-26,000)  =  -174,000.  With  the  oy  -  250,000  psi  material, 
yielding  will  not  occur  in  the  bore. 

To  f  .d  the  resiuuai  hoop  stress  upon  release  of  the  200.000  psi 
static  pressure  Pgoes  to  zero,  and  the  hoop  stress  is  ditr.'-r.iihed  by  the 
amount  =217,  OOOpsi.  Then 

Cfe  =  Oa  +  (4 

Cfe  =  133,000  +  (-110,000) 

Cfe  =  243,  000  psi 

Therefore  the  bore  diameter  at  rQ  =  a  =  3  inches  has  attained  a 
residual  hoop  stress  (compression)  of  243,000  psi  which  is  within  the 
required  250,000  psi  yield  strength  of  the  material.  The  pressure  vessel 
can  be  operated  safely  in  the  prescribed  limits  of  200,000  psi  static 
pressure  and  a  peak  shock  pulse  of  248,  280  psi  for  20  jx duration. 
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